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We review a striking array of recent experiments, and tlicir tlieoretical interpretations, on the 
superfiuid transition in ''He in the presence of a heat flux, Q. We define and evaluate a new set of 
critical point exponents. The statics and dynamics of the superfiuid-normal interface are discussed, 
with special attention to the role of gravity. If Q is in the same direction as gravity, a self-organized 
state can arise, in which the entire sample has a uniform reduced temperature, on either the normal 
or superfiuid side of the transition. Finally, we review recent theory and experiment regarding the 
heat capacity at constant Q. The excitement that surrounds this field arises from the fact that 
advanced thermometry and the future availability of a microgravity experimental platform aboard 
the International Space Station will soon open to experimental exploration decades of reduced 
temperature that were previously inaccessible. 
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Due to its extraordinary purity and insensitivity to ex- 
ternal perturbations, superfiuid ^He has long been the 
best system for accurate, detailed experimental inves- 
tigations of phase transitions and critical phenomena. 
The essence of superfluidity, however, lies in the dynam- 
ics of flowing helium, rather than in equilibrium proper- 
ties such as the specific heat and the superfiuid density. 
In recent years, considerable experimental attention has 
been focused on the behavior of liquid helium, near the 
lambda transition temperature, Ta, when a flux of heat, 
Q, is passed through it (Duncan et ai, 1988; Liu and 
Ahlers, 1996; Moeur et al, 1997; Day et at, 1998; Har- 
ter et ai, 2000). This article reviews the wealth of ex- 
citing new physical phenomena uncovered by these ex- 
periments and by the parallel theoretical investigations 
(Onuki, 1983; Haussmann and Dohm, 1992a; Haussmann 



and Dohm, 1994; Weichman et ai, 1998; Haussmann, 
1999a; Weichman and Miller, 2000). 

At sufficiently small Q, superfiuid ''He transports heat 
essentially without dissipation, by means of superfiuid- 
normal fiuid counterfiow. At higher Q, or as T ap- 
proaches Tx, this transport mode breaks down. In Sec. 
II we discuss the phase diagram of helium in the Q-T 
plane. In this plane, the superfiuid state is bounded by 
a transition curve, Qc{T), outside of which dissipative 
fiow takes over. We also introduce and evaluate a new 
set of critical exponents that arise as a consequence of 
superfiow. 

In Sec. HI we describe an inhomogeneous phase where 
the heat fiows through both normal and superfiuid re- 
gions, separated by an interface region that is neither 
normal nor superfiuid. In the normal region the heat 
flux produces a static temperature gradient. In the su- 
perfluid region, the heat flows at constant temperature. 
In the interface region, a transition between these types 
of behavior is mediated by fluctuations in a way that is 
not yet accessible either to experiment or to theory. 

A recurring theme throughout the article is the fun- 
damental limitation imposed by the Earth's acceleration 
due to gravity, ge, on the resolution of Earth-based ex- 
periments. Gravity produces a pressure gradient across 
the sample, leading to a variation in the local lambda 
transition temperature T\{z) with height z according to 



a,T,(z) = -7-, 7^1-27^. 

9e cm 



(1.1) 



For generality, the possibility of gravity g different from 
that on Earth is included. For g — ge the transition tem- 
perature therefore decreases by 1.27//K per centimeter of 
column height. If T is tuned so that, say, the center of 
the cell is at the local lambda point, the upper region will 
be superfluid while the lower region will be normal fluid. 
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The interface between, defined as the region over which 
the local properties differ significantly from those of a 
bulk system at the same local temperature and pressure, 
is about 0.1 mm wide on Earth (Muzikar and Giordano, 
1989). This figure also estimates the maximum possible 
critical correlation length in the system. 

The inhomogencity induced by ge causes the critical 
singularities to be rounded on a scale set by the height of 
the cell. Balancing gravity effects against finite size ef- 
fects, which enter if the cell is too small, the optimal cell 
height is about 0.3 mm, leading to rounding on a scale 
e = {T - Tx)/Tx ~ 10"'^. Present thermometry (Lipa 
et al., 1996; Day et ai, 1998) is capable of resolving re- 
duced temperatures e that are 3-4 orders of magnitude 
smaller than this. In essence there is a new frontier, con- 
sisting of decades of previously unavailable reduced tem- 
perature around Tx, that can only be explored in micro- 
gravity. Experiments measuring equilibrium specific heat 
(Lambda Point Experiment — LPE) (Lipa et al., 1996) 
and finite size effects (Confined Helium Experiment — 
CHEX) have recently flown aboard the space shuttle. 

Gravity also has a strong effect on the interface region 
discussed in Sec. IIL As we have seen above, the inter- 
face is compressed by gravity to a width of about 0.1 mm, 
too small to be stiidicd experimentally. The interface re- 
gion will be studied in an experiment (Critical Dynamics 
Experiment — DYNAMX) presently being prepared for a 
low temperature microgravity platform (Low Tempera- 
ture Microgravity Physics Facility — LTMPF) that is to 
be part of the International Space Station (ISS). 

Although gravity is detrimental to some measure- 
ments, there can also be interesting new physics when 
both gravity and heat current are present. In Sec. IV we 
describe the so-called self- organized critical (SOC) state, 
in which g and Q conspire to produce a new, essentially 
homogeneous, nonequilibrium state where the tempera- 
ture T{z) precisely parallels Tx{z) at a fixed Q-dcpendcnt 
distance. At larger Q, this state undergoes a transition 
from normal to superfiuid, with the temperature gradient 
supported by a stream of vortices in the latter. 

Finally, in Sec. V we discuss the specific heats under 
heat flow, predicted (Haussmann and Dohm, 1994; Chui 
et ai, 1996; Haussmann and Dohm, 1996), and subse- 
quently confirmed (Harter et al, 2000), to be enhanced 
above the equilibrium form, with a singularity predicted 
at the phase boundary Qc{T). This experiment too ulti- 
mately needs to be performed at very low Q, and the re- 
quired microgravity version is currently being proposed. 



II. NONEQUILIBRIUM CRITICAL PHENOMENA AND 
SCALING 

We begin by presenting the basic mathematical back- 
ground within which the physical phenomena we discuss 
in later sections is most clearly described and understood. 



The language of phase transitions and critical phenom- 
ena, under both equilibrium and nonequilibrium condi- 
tions, is that of critical exponents, scaling relations, and 
scaling functions, which quantify the scale invariant na- 
ture of a system near its critical point. In order to keep 
the discussion at as elementary a level as possible, we will 
introduce the theory in close analogy to that appropriate 
to equilibrium classical magnetic and liquid-vapor criti- 
cal points. The near-dissipationless nature of superflow 
in '*Hc makes this analogy especially fruitful in describing 
nonequilibrium superfluidity. 



A. Superfiuid counterflow 

Under conditions where a uniform heat flux Q is ap- 
plied to superfluid ^He, a thermal counterflow is created 
in which the normal fluid moves along with Q at veloc- 
ity u„, and the superfluid moves oppositely with velocity 
Us. For sufliciently small heat current not too close to Tx 
[more precisely, one rcqiures Q < Qc{T), with Qc{T) de- 
flned in Fig. 2 and (2.20) below] the temperature T <Tx 
is essentially uniform and the heat flux is 



Q = TSpUn, 



(2.1) 



where S is the entropy per unit mass, p = Ps + Pn is the 
total mass density, composed of superfluid and normal 
fluid parts, and u„ is the normal fluid velocity. Since 
there is no net mass flow, one has j„ = — js, where 
js = PsUs and j„ = p„u„ are, respectively, the superfluid 
and normal fluid mass current densities, and is the 
superfluid velocity. Close to Tx one finds experimentally 
S ^ Sx = 1.58J/gK, p, « po\e\i, with po ~ O.STg/cm^, 
critical exponent ( ^ 0.671, and pn — p — 0.14g/cm^. 
Thus, in experimentally motivated units, 



Us « 7.9 X 10"^ 
Un w 2.1 X 10-'^ 



Q 



10 



-6\ C 



l/xW/cm2 \ |e 

Q 

lAtW/cm2 



cm/s (2.2) 
cm/s. (2.3) 



At experimentally accessible values Q = l/xW/cm^ and 
|e| = 10~^, Us is a modest 80 /xm/s and u„ is nearly four 
orders of magnitude smaller. 



B. Thermodynamic formalism 

The isothermal condition allows an effective thermody- 
namic description of the finite Q state (Hohenberg and 
Martin, 1965). Although nonequilibrium scaling does not 
rely on this, a more intuitive description of a number of 
phenomena special to superfiuidity results, and so it is 
worthwhile presenting the theory in this context. 
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In the frame of reference moving with the normal fluid 
(in which there is no heat flow), the differential of the 
free energy density F(T, Ug) may be written, 



(a) h 



dF = -SdT + J, ■ dVs 



in which U, = u, — u„ and 



dF 



(2.4) 



(2.5) 



with the second equality serving as a definition of 
Ps{T, Us) when Us is not small. Near T\ the smallness of 
u„ implies that Us ~ Us and Jg — js- In this frame one 
has the defining relation = (7i/m)V(?i, where m is the 
*He atomic mass, and 



777 

0(r) = -U,.r 



(2.6) 



i4> 



is the phase of the superfluid order parameter tp = lifjle 
Thus, V' rotates in the complex plane as one moves along 
the direction U^, generating a kind of helical structure. 
The effective thermodynamic description (2.4) relies on 
the existence of a time-independent tp in the presence 
of a finite phase gradient. In fact, thermally nucleated 
phase slips in the helical structure (interpreted as tunnel- 
ing between different metastable local minima of the free 
energy) lead to small temperature gradients and decay 
of superflow. However, the decay time is extremely large 
at low heat currents and temperatures not too close to 
Tx (Langer and Fisher, 1967). The thermodynamic de- 
scription (2.4) is valid on time scales smaller than this. 
As Us increases, both the order parameter magnitude 
ItpiUsjT)] and the superfluid density are suppressed rel- 
ative to their equilibrium values at J7s = 0, Q = 0. 

Variations at constant Q = —TSJs are most conve- 
niently performed by defining the Legendre transformed 
free energy $(T, Js) = F — Jg - Ug with differential 



= -SdT - Us • dJs 



(2.7) 



Close to T\, where TS ~ TxSx is essentially constant, 
variations at constant Q are asymptotically the same as 
those at constant Js. For example, the specific heat at 
fixed Q may be taken as 



Cq 



dS_\ 



(2.8) 



where S{T,Js) = -{d^/dT)j^. With the above obser- 
vation in mind, we shall henceforth treat and Q as 
differing only by a constant factor. 



C. Nonequilibrium scaling 

The fact that Us and Js (or Q) may be treated as ther- 
modynamically conjugate variables has important conse- 
quences for the structure of the thermodynamic functions 



1st order 
coex. line 



(b) Q 



dissipative phase 
..,Qc(T) 





FIG. 1. Compaxison of phase diagrams in (a) conventional 
T-h and T-m spaces and (b) superfluid T-Q and T-Us spaces. 
Here, mo{T) is the spontaneous magnetization, or the liquid 
density minus the critical density. The critical lines QdT) 
and Us,c{T), enclose the region of stable superflow, which 
corresponds also to the region of validity of the thermody- 
namic description, and like the thermodynamic description 
itself are sharp only in the absence of phase slips. The dif- 
ferent shapes of those curves near T\ arc determined by the 
fact that Aq > 1, while /3q < 1 [see equations (2.18)-(2.21) 
below]. The nature of the critical behavior as Q approaches 
Qc from below, and the nature of the inhomogeneous dissi- 
pative phase for Q > Qc, will be discussed in later sections. 
A major diff'crcnco between (a) and (b) is the lack of a first 
order line below in the latter. Thus, h vanishes through- 
out the phase coexistence region in the T-m plane, whereas 
Q oc Js = PsUs varies continuously throughout the superfluid 
phase in the T-Us plane. 



near Tx (Fisher, 1973). We proceed by analogy with 

the conjugate variables h and m at conventional criti- 
cal points, where h is the external magnetic field and 
m the magnetization at a Curie point, or h is the dif- 
ference from critical pressure and m the difference from 
critical density at a liquid- vapor critical point. There 
the free energy A{T, h), analogous to has differential 
dA = —SdT — mdh and its singular part As obeys an 
asymptotic scaling form [see, e.g., (Fisher, 1983)], 



As(T,/i)=^o|e|'""^ 



Dc^h 



(2.9) 



in which a is the specific heat exponent, A is the "gap 
exponent," A{x) is a universal scaling function, and Eq, 
Do are nonunivcrsal scale factors, specified uniquely via, 
say, the normalizations ^(0) = ^'(0) = 1. There are 
actually two scaling functions, A±{x) for ±e > 0, but we 
shall primarily be interested in the ordered phase e < 
and consider only A- = A. 

Similarly, the singular part of the superfluid free en- 
ergy, $s, is expected to obey the scaling form, 
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^siT,Q)=Ao\ef-''Y 



(2.10) 



in which a ~ —0.013 (Lipa et ai, 1996) is again the 
usual equihbrium specific heat exponent, Aq is the gap 
exponent for Q, Y{y) is the e < universal scaling func- 
tion, and Aq, Qq are nonuniversal scale factors, specified 
uniquely via, say, the normalizations 5^(0) = Y"{0) = 1 
[Y{y) must be an even function of y due to the obvious 
symmetry under sign reversal of Q] . 

The derivative of (2.9) with respect to h yields 



m{T,h) = -EoDo\efA' 



Doh 



(2.11) 



where the prime denotes derivative with respect to ar- 
gument and the order parameter exponent /3 obeys the 
scaling relation (3 = 2 — a — A. The second h derivative 
yields the order parameter susceptibility (compressibility, 
in the case of a liquid-vapor critical point) 



xiT,h) 



(2.12) 



with 7 = q:+2A — 2, which then yields the famous Essam- 
Fisher scaling law a + 2/3 + 7 = 2 (Fisher, 1983). 

Similarly, the derivative of with respect to Jg yields 
Us in the form (Haussmann and Dohm, 1992b) 



where Bq = AqTxSx/Qq, and one has the generalized 
order parameter exponent scaling relation 



/3Q = 2-a-AQ. 



(2.14) 



The equilibrium superfluid density enters the free energy 
F via a term AFs = \psUa for small Us- The Legendre 
transform yields a term A$s = —Jljlps for small Jg, 
and the inverse of the superfluid density now appears in 
the theory as a generalized susceptibility: 



1 



P«(r,g = o) 



T,J,=0 



dUs 
dJs 



T,Js=0 



RqY" jO) 
|e|T« 



(2.15) 



where Rq = Ao{T\S\/Qo)^ and the generalized sus- 
ceptibility exponent jq obeys the Essam-Fisher rela- 
tion (Fisher, 1983), 



7Q = C = 2Aq 
a + 2/3q + 7Q 



a - 2 



(2.16) 
(2.17) 



Generally the gap exponent is independent of a and must 
be separately determined. However, in the superfiuid 
problem the Josephson relation [see, e.g., (Pfeuty et al., 



1974) and references therein] yields C = 2 — a — 2j/ = 
{d — 2)1/. Here ^ ~ ^o/leT with, in dimension d = 2>, u = 
C ^ 0.671 and ^0 3.4A (Goldner and Ahlers, 1992), de- 
scribes the divergence of the superfluid coherence length, 
^ = {w?kBT /H? psY^^'^~^\ and the second equality fol- 



lows from the hyperscaling relation 2 — a 
1973; Fisher, 1983). We therefore identify 



2-a-v={d- l)v. 



dv (Fisher, 



(2.18) 



This scaling law implies that Q scales with the cross- 
sectional area S^"^"^ of a correlation volume Q becomes 
significant when the power incident on a correlation area 
is of order Qo^o~^* Prom (2.14) one obtains 



(2.19) 



This relation has the interpretation that {m/h)Us, which 
has dimensions of inverse length, scales with the inverse 

correlation length the phase gradient has a signifi- 
cant effect when its wavelength 2Trh/m.Us becomes com- 
parable to ^. Note that, more typically, one begins with 
the latter assumption and reverses the above argument 
to derive the Josephson relation. 

Just as arbitrarily small h smears the singular behavior 
near a Curie point, one expects an arbitrarily small Q to 
either smear or drastically alter the lambda point critical 
behavior. This is consistent with Aq > 0, implying that 
the scaling argument y diverges as [f] at any finite 
Q, serving to define Q a,s a relevant perturbation to the 
lambda point (Fisher, 1973). In Fig. 1 a schematic of 
the expected phase diagram is shown, contrasting it with 
that for a conventional critical point. The lines Qc{T) 
and Us,c{T) are the boundaries beyond which superflu- 
idity breaks down: for Q > Qc, due to suppression of 
the order parameter and superfiuid density, the heat cur- 
rent is too large for the superfluid to support isothermal 
heat transport, and a nonequilibrium phase transition to 
a new dissipative phase occurs. The nature of this phase 
will be discussed in Sec. Ill below. In Fig. 2 we sketch the 
isothermal equations of state for the conventional and su- 
perfluid systems. At a fixed subcritical temperature the 



*One may estimate Qo as follows. Two scale factor uni- 
versality yields a form Cs = ksiR^/^)'^ /a\£\'^ for the sin- 
gular part of the equilibrium specific heat below T\, where 
the hyporunivcrsal ratio ~ 0.90 in d = 3 (Hohenberg 
et al, 1976). This must match the scaling form (2.10) at 
Q = and determines Ao|e|^"" = [kBTx/a{l - a)(2 - 
a)]{R^/^)'^ [with the choice y(0) = 1]. Equating the 
quadratic terms Ao\e\'^'°' {Q/Qo\e\'^'^f = -Js/Ps [with the 
choice F"(0) = 1], one obtains finally Qo = T\S\\/—Aopo, 
and Qaif~^ = mksT^SxR'^/^ /hy/ -a(l - a){2 - a) (using 
po = m?kBTx/h^£,Q~^). Substitution of experimental num- 
bers yields Ao ~ -21J/cm^ and Qo ^ 30kW/cm^. 
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FIG. 2. Comparison of the isothermal equation of state (a) 
in the conventional h-m plane and (b) in the superfluid Q- 
Us plane. At /i = there is a first order transition between 
up and down magnetized states, or between liquid and va- 
por states. As increases, \rn{T,h)\ increases as well. As 
IQI increases \Us\ increases until the superfluid breaJts down 
at IQI = Qc(T), at which point Us{Q,T) has a square root 
cusp. At this same point the specific heat at constant Us also 
has a square root cusp (Haussmann and Dohm, 1994), while 
the specific heat at constant Q has an inverse square root 
divergence (Chui et al., 1996). 

conventional system displays the usual first order jump 
in m at ft, = between up and down magnetized states, 
or between liquid and vapor states. The superfluid sys- 
tem displays a continuous variation of Us with Q, but 
at IQI = Qc{T) encounters the boundary between super- 
fluid and dissipative phases. At this point Us{Q) has a 
square root cusp (Haussmann and Dohm, 1994), corre- 
sponding to a strong suppression of ps and a divergent 
susceptibility {dUs/dJs)T- This singularity is exhibited 
in the scaling form (2.10) as square root cusps in Y{y) at 
some finite value y = dszyc- This yields the predictions 

Qc{T) « QoyM^'' (2.20) 
Us AT) « BoY'{y,)\ef'i « Qc(T)i/^«, (2.21) 

in which 

Sq=Aq/Pq = 1 + jq/Pq- (2.22) 

The latter relation generalizes the Widom scaling law 
m{Tc,h) oc h^^^ with d = A//3 = 1 + 7//? which fol- 
lows from (2.11) in the limit e ^ [the asymptotic be- 
havior A'{x) ~ x^/^ as a; ^ 00 is required for consis- 
tency] (Fisher, 1983). From (2.16) and (2.19) one obtains 
explicitly Sq = d — 1. As will be discussed in Sec. V, the 
specific heats at constant Us and Q are also singular at 
Us,c{T) and Qc{T) (Haussmann and Dohm, 1994; Chui 
et a/., 1996; Haussmann and Dohm, 1996). 
The coherence length itself obeys a scaling form 

«r.Q) = a|er'H(^), (2.23) 

with S(0) = 1. Considering variations of ^ at a fixed 
value of y = Q/Qol^]^'^, one may write 

^ =Co(Q/Qo)"''«y'^«S(y) (2.24) 
UQ = u/AQ = l/{d-l), (2.25) 



exponent general d d = 3 

V' 1/2 < i^{d) < 00 0.671 

a 2-dv -0.013 

C 2- a-2v = {d-2)v 0.671 

Aq 2-a-u={d-l)u 1.342 

(3q V 0.671 

7Q C 0.671 

l+7Q//?Q=rf-l 2 

VQ i//AQ = l/(d-l) 1/2 

TABLE I. Equilbrium and nonequilibrium critical exponents 
and their values. Note that no exact expression for the corre- 
lation length exponent v{d) is known, aside from the bound- 
ary values ;^(2) = 00 and v(d > 4) = i, but all other expo- 
nents are either given exactly or in terms of v. There are other 
independent exponents (e.g., the critical correlation exponent 
T] ~ 0.02), but they happen not to be involved in any of the 
experiments we discuss. Expressions in the center column 
involving d explicitly require hyperscaling, and are therefore 
valid only for d < 4. For d > 4 the exponents all stick at 
their mean field values, and may be determined by substitut- 
ing V = i, a = into the non-explicitly d-dependent forms 
of the scaling relations. 

serving to define an exponent describing the characteris- 
tic variation of ^ with Q. Attempts to verify the scaling 
relations (2.20) and (2.25), with Ag = 2v and i/q = 1/2 
in = 3, will be discussed in later sections. 

In Table I we summarize the various critical exponents 
we have defined, along with their values in general di- 
mension d and in d = 3. 

As a final comment, we note that verifications of cer- 
tain dynamical scaling laws in ^He are complicated by 
a fundamental problem associated with the existence of 
very slowly convergent, nonasymptotic dynamic effects. 
Scaling forms like (2.10) implicitly assume that |e| is suf- 
ficiently small that all irrelevant scaling variables may 
be ignored. It turns out that a certain Wegner expo- 
nent ui controls the relevance of mass and heat diffu- 
sion on the dynamic critical behavior (Hohenbcrg and 
Halperin, 1977) through an additional scaling variable 
w = 7o|e|'^, where 70 is a parameter of order unity in 
the Model F equations (Hohenberg and Halperin, 1977) 
(describing the near-critical dynamics of ^He) which cou- 
ples entropy and density fluctuations to order param- 
eter fluctuations. The variable w does not appear in 



'The Hohenberg-Halpcrin Models A-J classify the different 
dynamical universality classes encountered most commonly in 
physical systems. The dynamical behavior depends qualita- 
tively, for example, on whether or not the order parameter is 
a conserved density, and whether or not it is coupled to other 
conserved densities (like mass or energy) . The Model F equar- 
tions axe the simplest possible description of a non-conserved 
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any equilibrium scaling function (where such fluctua- 
tions may be completely "integrated out" of the partition 
function), but does appear in those involving dissipative 
transport, e.g., that of the thermal conductivity (Hauss- 
mann and Dohm, 1992a). Since io ~ 0.008 (Dohm, 1991) 
is very small at the superfluid transition this variable van- 
ishes extremely slowly as |e| ^ 0: e = 10^^^'^ leads only 
to lel*^ = 0.1. The scaling function then has a very slow 
parametric dependence on w, leading to "quasi-scaling" 
(Haussmann and Dohm, 1992a) of the heat conductivity, 
and hence to an apparent slow variation of the associated 
critical exponent with w. Measured dynamic critical ex- 
ponents affected in this way may typically be expected 
to differ from their true asymptotic values by 10-20% 
(Duncan et al., 1988). 

The scaling phenomena considered here, however, are 
mainly those associated with properties of isothermal su- 
perfluids, which though out of equilibrium, are never- 
theless argued to behave as equilibrium thermodynamic 
systems. To the extent that this is true (i.e., to the ex- 
tent that vortex excitations can be ignored), the Model 
F equations again produce a thermodynamic-type par- 
tition function with entropy/density fluctuations inte- 
grated out, and finite Q enforced by an imposed uniform 
helical twist in the order parameter (Haussmann and 
Dohm, 1994). The scaling variable w will again not ap- 
pear, and one expects that quasi-scaling will be absent in 
(2.10) and from all quantities derived from it, and hence 
that the exponents in Table I will exhibit experimentally 
their predicted values. On the other hand, the phase 
boundaries in Fig. 1(b) are defined by the onset of dissi- 
pative transport, with divergent fluctuations in the local 
heat current as they arc approached. Such fluctuations 
lead to vortex creation, decay of supcrflow, breakdown of 
the effective equilibrium description, and the reappear- 
ance of the variable w. In some sense w must control the 
"fuzziness" of the boundary Qc(T), and an experimen- 
tal test of the relations (2.20) may exhibit quasi-scaling 
even if, for y sufficiently less than j/c, the scaling function 
Y{y) does not. This issue is potentially testable by the 
specific heat measiircments discussed in Sec. V. Unfor- 
tunately, as exhibited in Fig. 9 below, present data are 
constrained to lie sufficiently far below Qc{T) that both 
asymptotic and nonasymptotic forms fit equally well. 



two-component order parameter (the real and imaginary parts 
of V) coupled to heat and mass flow, and are believed to 
provide an asymptotically exact description of the universal 
critical dynamics of superfluid '*He. 



III. THE NONEQUILIBRIUM SUPERFLUID-NORMAL 
INTERFACE 

A. Interface statics 

Consider a cylindrical cell with a heat current Q driven 
along its axis, labeled by coordinate z, in which the up- 
stream endwall, 2 = 0, has T(0) > T^. In the normal 
phase heat is transported by thermal conduction, which 
at sufficiently small Q is described by the Fourier law 

Q = -Kd,T, (3.1) 

where k{T) « kqc"'' is the thermal conductivity, pre- 
dicted to diverge at Tx, consistent with its infinite value 
at all T < T\. Experimentally one finds (Dingus, et 
al, 1986; Tam and Ahlers, 1985; Moeur et al, 1997) 
Ko — 12/xW/cm^ and /x ~ 0.44. In a very tall cell in 
which T varies substantially along its length, one may 
view (3.1) as locally valid with k(z) = K[r(z)] so long 
as T[z) remains sufficiently far above Tx. "Sufficiently 
far above Ta" may be quantified by the condition that 
the temperature drop across a coherence length be much 
smaller than the deviation from Tx: 

^|VT|«T-Ta ^ (3.2) 

Putting in ^He parameters, this requires 

e>6xl0-^f— — , r=- ^ ^0.81 

Vl/xW/cm^y l + v-n. 

(3.3) 

As z increases, T{z) will eventually violate (3.3), and one 
enters a region of nonlinear heat transport. In effect, k 
becomes a strong function of Q in this region. Moreover, 
as illustrated in Fig. 3, at some position zq, T{zo) = Tx, 
and the system enters the superfluid phase iox z> z^: a 
nonequilibrium superfluid-normal fiuid interface is gen- 
erated. Far downstream from this interface, T[z) levels 
out at a temperature Tao{Q), one of the nonequilibrium 
thermodynamic states discussed in Sec. II above. Corre- 
spondingly, the order parameter magnitude |V'(-2)|, which 
effectively vanishes for z < zq, grows in the interface re- 
gion and saturates at a value |V'oo(Q)| for z':^ z^.^ 

The interface is the region over which the mode of heat 
transport converts from conduction to superfluid coun- 
terflow. Its characteristic width ^(Q), through which T 



*In fact T{z) continues to decrease slowly with z in the 
superfluid phase due to rare phase slip events (Haussmann, 

1999a; Haussmann, 1999b), so Too{Q) and I'tpooiQ)] arc not 
sharply defined, but we shall not dwell on this complication 
here. 
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FIG. 3. Scaled temperature and order parameter profiles 
M(Z) and ^(^), computed within the mean field approxima- 
tion (Onuki, 1983; Weichman et al, 1998). Here Z = z/l{Q), 
1*1 oc WiQ) and M oc (T - To)f{Q), where To is the mean 
field transition temperature and 1{Q) oc Q~^^^ is essentially 
the correlation length (2.24) in the mean field approximation. 
With these scalings the profiles are independent of Q. 



and vary substantially (Fig. 3), is expected to scale 
according to (2.24). We may estimate this width semi- 
quantitatively by defining a reduced temperature scale 
eg through equality in (3.2) and (3.3), and defining 



Q 



1 fiW/c 



/im, riy ~ 0.54. 

(3.4) 



The scaling relation (2.25) predicts rv = vq (= 1/2 in 
d ~ i), requiring the scaling relation = 1 — C (~ 1/3 in 
d = 3) (Hohenberg and Halperin, 1977). The difference 
between experimental and theoretical values is presum- 
ably due to the quasi-scaling effect described at the end of 
Sec. II. Quantitative renormalization group based predic- 
tions for the full temperature profile T{z,Q), explicitly 
including the quasi-scaling phenomenon (Haussmann and 
Dohm, 1991; Haussmann and Dohm, 1992a), are limited 
to the normal fluid region z < zq where the vanishing 
of "0 greatly simplifies the calculations. Less quantita- 
tive predictions for the entire profile are limited to mean 
field (Onuki, 1983) (which ignores critical fluctuations) or 
large- (Haussmann, 1999a; Haussmann, 1999b) (which 
replaces the single complex order parameter with an A^- 
component vector) approximations. 

An experimental measurement of the temperature pro- 
file r(z, Q) within the interface region for a sequence of 
different Q would allow a detailed exploration of near- 
critical, nonlinear heat transport. It would also pro- 
vide an experimental test of the scaling prediction (2.25). 
Such a measurement requires a regime in which ^{Q) is 
larger than the width W of the thermometer. Present 



technology places a limit W > 50 /J,m. From (3.4), heat 
currents below 1 /iW/cm^ are required. Present tech- 
nology allows controlled heat currents down to about 
1 nW/cm^, leading to ^ w 1 mm, so at first sight such 
an experiment appears feasible. 

In fact. Earth's gravity ge places a fundamental limit 
on the maximum possible interface width. As described 
in Sec. I, gravity produces an equilibrium [Q — 0) 
superfluid-normal fluid interface with width of order 
e(.9e) = 100 From (3.4), f(Q) = ^{g,) for Q ^ 

Qg = 70 nW/cm^. For Q of order Qg, gravity will begin 
to have a strong effect on the nonequilibrium interface 
(Day et ai, 1998), and for Q <^ Qg the heat current 
will be a small perturbation on the equilibrium interface. 
Thus ^(Q) will saturate at £,{ge) as Q ^ and the regime 
^ W is unattainable on Earth. 

For this reason the critical dynamics experiment (DY- 
NAMX) is currently being prepared for the micro- 
gravity environment of the International Space Station, 
with flight planned for 2004. Heat currents as low as 
5 nW/cm^ will be used. The temperature profile T(z, Q) 
will be measured to subnanoKelvin resolution by slowly 
moving the interface past a fixed thermometer. The lat- 
ter is accomplished by removing heat from the down- 
stream end of the cell slightly more slowly than it enters 
the upstream end, with the effect that the normal phase 
slowly invades the cell, thus translating the interface. 



B. Transition from thermodynamic to interface state 

We have discussed three possible classes of steady 
nonequilibrium states: (a) conducting normal states 
with static temperature profile determined by (3.1); (b) 
isothermal thermodynamic states with steady superfluid 
counterflow determined by (2.1); and (c) states with a 
nonequilibrium interface forming a "conversion bound- 
ary" between states of type (a) and (b) . Transitions be- 
tween (a) and (c) occur continuously: if state (a) is cooled 
to the point where the downstream endwall temperature 
passes through T\, an interface will form out of that end- 
wall and steadily move upstream. Conversely, as heat is 
added to state (c), the interface will move downstream 
until the interface disappears into the endwall. 

Similarly, the transition from (c) to (b) is continuous. 
The interface will disappear into the upstream endwall as 
heat is extracted, yielding an isothermal thermodynamic 
state at temperature Too{Q)- Further extraction of heat 
will cause the temperature to drop below Taa{Q). 

The nature of the transition from (b) to (c) is less clear. 
At issue is whether the bulk superfluid, in the absence of 
an interface, recognizes Too{Q) as a special temperature. 
Within the mean field approximation, the answer is no 
(Onuki, 1983): the thermodynamic state is stable up to 
a temperature Tc(Q), with Tx > Tc{Q) > Toa{Q), whose 



7 



M 

-1 .6 



-2.6 





\ 

'10 




G=0 


9 

^ 

\ 

\ 


\ 

\ 

I 

\ 






\ 

^ \ 

\ \ 
\ \ 


\ 

\ 

\ 






\ \ 
-\ \ 
7 \ \ 

-\\\ 


\ 

\ 

\ 

\ 





- 




^ 














































■r---- 

















10 15 20 — 25 30 



FIG. 4. For bulk scaled temperature M > Mc (see caption 
to Fig. 3) the uniform superfluid state becomes unstable to a 
state with an interface (Weichman and Miller, 2000). Curves 
1-6 (solid) show consecutive stable (time independent) super- 
fluid temperature profiles obtained in the mean field approxi- 
mation upon slowly varying the right-hand wall temperature. 
Curves 7-10 show snapshots of the time evolution, computed 
within a simple one-dimensional model, initiated by a series 
of phase slips (not shown) after the right-hand wall tempera- 
ture is raised slightly above that of curve 6. The key feature 
is the net temperature drop (from Mc to Moo) observed in the 
bulk superfluid. Figure reprinted from (Weichman and Miller, 
2000) courtesy of the Journal of Low Temperature Physics. 

functional inverse Qc{T) we might identify with the 
boundary in Fig. 1(b), at which the superfluid density is 
suppressed to the point where it is incapable of support- 
ing the heat current. As shown in Fig. 4, when the ther- 
modynamic state is heated above TdQ) a complicated 
dynamics results, with the system finally settling down 
into a state with an interface. Since Too{Q) < Tc{Q), the 
superfluid side actually cools, with a compensatory heat- 
ing of the normal side. The final position of the interface 
is determined from energy conservation. 

The question of whether or not a first order transi- 
tion from (b) to (c) survives in some form beyond the 
mean field approximation is subtle, and is probably a 
question of time scale. The same thermal micleation of 
vortices that leads to a small temperature gradient in the 
superfluid phase presumably also leads to a continuous 
(b) to (c) transition, with the interface entering contin- 
uously from the upstream endwall, if the experimental 
heating rate is infinitesimally slow. Any finite heating 
rate may, however, allow a "superheating" of (b), induc- 
ing an apparent first order transition to (c) nucleated by 
a random vortex creation event. Such an effect would 
be analogous to superheating and supercooling effects at 
conventional first order transitions, the boundary Qc{T) 
being analogous to a spinodal line (defined roughly as 
the temperature at which a local free energy barrier sep- 
arating the superheated or supercooled metastable state 



from the true equlibrium state disappears, allowing rapid 
conversion to the thermodynamically stable phase). 

Experimental investigation of the (b) to (c) transition 
is complicated by boundary effects. The singular Kapitza 
resistance leads to an additional heating of the upstream 
endwall, which then acts as a nucleation point for vortices 
even when the temperature of the bulk superfluid is still 
below Too(Q) (Harter et al., 2000). Clever cell designs 
that reduce the heat current near the endwall below that 
in the bulk may eventually allow experimental observa- 
tion of a superheating effect. 

C. Interface dynamics 

Despite the greatly reduced gravity, the Space Sta- 
tion environment is less than ideal for other reasons. 
Vibrational noise ("g-gitter") exists at a level of about 
10~^(7e, and one might worry that such noise might cou- 
ple strongly to the interface, and perhaps destabilize it. 
Understanding the effects of acceleration noise requires 
an understanding of (a) the free dynamics of the inter- 
face when it is perturbed away from its steady state, e.g., 
whether or not the interface is even dynamically stable, 
or whether small perturbations might undergo some kind 
of "dendritic growth", and (b) the manner in which vi- 
bration, or other perturbations, couple to, and perhaps 
amplify this motion (Weichman et al., 1998). 

Suppose that a slow variation, zo = zo{r), with the 
transverse coordinates r = {x, y), is applied to the inter- 
face position and is subsequently released. The problem 
is to derive an equation of motion for zo{r,t). In order 
to motivate the interesting physical questions, it is useful 
to consider first the analogous problem of the motion of 
an equilibrium interface between up and down domains of 
an Ising fcrromagnet, with diffusive "Model A" dynamics 
(Hohenberg and Halperin, 1977). Surface tension acts as 
a restoring force against perturbations away from a flat 
interface, and one may derive a diffusion equation for the 
relaxation of long wavelength perturbations of the inter- 
face position: 

(5t-DV2)^o(r,t) =??(r,t), (3.5) 

in which the value of D depends on the microscopic pa- 
rameters in the model, and the driving term r/, vanish- 
ing for free relaxation, has been included for complete- 
ness. Thermal fluctuations in the spins lead to a Gaussian 
white noise form for r] with correlator {rj{r,t)r}{r' ,t')) = 
Co<5(r — r')6{t — t') where C^q depends on microscopic pa- 
rameters and on temperature. Using this form one may 
compute the equal time variance in the interface position 

C(r-r') = (N(r,t)-^o(r',t)]2) 

= (3.6) 
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where 0,0 is an atomic length scale. One sees that fluctu- 
ations in the interface position diverge logarithmically 
with separation: a famous result known as interface 
roughness, encountered most frequently in discussions of 
crystal facet shapes. The interface is locally stable and 
flat, but globally wanders arbitrarily large distances. 

The result (3.6) shows that there are subtle physical 
issues, lying beyond the much simpler question of stabil- 
ity, arising because the interface breaks the continuous 
translation invariance of the system. Whenever a con- 
tinuous symmetry is broken, a Goldstone mode (namely, 
a slow dynamical mode arising when the local value of 
the broken symmetry varies slowly in space) is gener- 
ated, corresponding in this case to very slow relaxation 
of long wavelength perturbations of the interface posi- 
tion, zo oc e»k r-A(k)t^ ^-^^i^ (3 5) yielding A = Dk'^. Such 
long wavelength "modes" are highly susceptible to ther- 
mal or externally generated noise spectrum C{k), and it 
is the convergence at small k of the integral 

that determines whether or not the interface is rough. 
For the case ({k) = (o, (3.7) diverges logarithmically, 
which same divergence is reflected in the correlator (3.6). 

The physics of the nonequilibrium superfluid-normal 
interface is very different from that of the equilibrium 
magnetic interface. As described, the dynamics of the 
latter is purely diffusive, with the two bulk phases on 
either side of it containing no slow dynamical modes of 
their own (since the order parameter is non-conserved 
and is not coupled to any other conserved field). In 
contrast, transport on the normal side of the ^He inter- 
face is controlled by slow heat diffusion, and transport 
on the superfluid side, controlled by superfluid counter- 
flow, is essentially ballistic, and the dynamics of the in- 
terface itself therefore involves a very intricate coupling 
of three different slow dynamical modes: normal phase 
diffusion, superfluid phase counterflow, and the broken 
translational symmetry Goldstone mode of the interface 
itself. To leading order one finds an equation of motion 
(Weichman et ai, 1998) 

{d^-c^V^)zo{r,t)=r]{r,t), (3.8) 

so that the interface supports traveling wave excitations, 
with a well defined speed c{Q) of the same order as 
the bulk second sound speed (the speed with which per- 
turbations in the order parameter travel on the supcr- 
fiuid side of the interface). At next- to- leading order one 
finds that these excitations are singularly damped, with 
A = ick+Dk^/"^ , so that RcA oc k'^/'^ rather than fc^ (as for 
bulk second sound waves) at small k. Positivity of Re£) 
establishes intrinsic dynamical stability of the interface. 
One may interpret the enhanced damping as arising from 



the waves on the interface "rubbing" up against the nor- 
mal phase. Moreover, one finds that thermal fiuctuations 
enter via a spectrum (^{k) oc A;^/^ for r/, vanishing strongly 
as fc ^ 0. The physics of this result is related to the fact 
that interface motion is a cooperative phenomenon, in- 
volving evanescent dynamics of the superfluid order pa- 
rameter to a depth ~ fc""^/^ scaling as the 3/2 power of 
the wavelength.^ The microscopic thermal noise, which 
is white, must then be averaged over a similar volume 
to obtain its net effect on the mode, leading eventually 
to the greatly reduced (^{k) above. In contrast, the Ising 
interface moves by local spin flips and the microscopic 
noise is averaged only over a microscopic region of width 
oq and C remains white. 

The net result of the analysis is that the integral (3.7) is 
strongly convergent at fc = 0, and the superfluid-normal 
interface is globally flat. This is good news for the Crit- 
ical Dynamics Experiment (DYNAMX), where a rough 
interface would have led to substantial smearing of the 
temperature profile on a scale varying as the logarithm 
of the cell cross-section. 

Vibrational acceleration noise couples to the interface 
in the same way that Earth's gravity does, through the 
variation in the local Tx with pressure. A slight change in 
T\ will cause the interface to translate (for acceleration 
normal to the interface) or tilt (for acceleration parallel 
to the interface) slightly, but if the change is oscillatory 
it coiild resonate with one of the intcrfacial second sound 
modes, leading to a rapid growth in the interface mo- 
tion. A detailed examination of the expected frequency 
spectrum of the Space Station g-gitter, together with the 
discrete spectrum of standing wave modes allowed in the 
experimental cell within the planned temperature and 
heat current range, shows that such resonances may in- 
deed occur, but that the singular damping is sufficiently 
strong that the interface oscillation amplitude should sat- 
urate at acceptably low levels (Chui et ai, 1997). 

The existence of the interfacial second sound mode has 
yet to be tested experimentally. This might bo accom- 
plished by applying a sequence of heat pulses to the cell 
sidewall near the interface and detecting a response at 
the opposite sidewall. Scattering of bulk second sound 
pulses off the interface, with detection of the reflected 
pulses, might also provide interesting information about 
the coupling of bulk and interface modes. 



^The evanscent depth on the normal side is much smaller, 
scaling as For comparison, gravity waves on fluid sur- 

faces yield fluid motion to a depth ~ /c~^ proportional to the 
wavelength. 
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IV. THE SELF-ORGANIZED CRITICAL STATE 

We have so far discussed phenomena in which optimal 
conditions occur in the absence of gravity. It transpires 
that there is a very interesting phenomenon in which 
gravity and heat current combine to produce a new type 
of dynamical state. The so-called self-organized critical 
(SOC) state occurs in a cell which is heated from above, 
so that g and Q are parallel.** 

Gravity depresses the lambda point T\(z) with in- 
creasing depth according to (1.1), while heat current 
leads to decreasing temperature T{z) with depth accord- 
ing to (3.1). The reduced temperature e(z) = [T{z) — 
T\{z)]/T\^Q, where T\fi is, say, the bulk transition tem- 
perature if gravity were absent (and hence approximately 
the transition temperature at the top of the cell), con- 
tains a competition between these two effects. One might, 
in fact, imagine tuning Q in such a way that e(z) is inde- 
pendent of z (Onuki, 1987): the sample would apparently 
exist in an essentially homogeneous near-critical state. In 
fact, it was argued (Machta et ai, 1993; Ahlers and Liu, 
1996) that the system actually "self-organizes" T{z) in 
order to enforce a uniform e. Assuming the validity of 
the Fourier law (3.1), e must be uniquely defined by 

^(esoc) ^Q/I^.TaI, (4.1) 

where esoCj ^ function of the ratio Q/g, is the reduced 
temperature of the new state. This state has recently 
been observed experimentally (Moeur et ai, 1997). 

As an aside, we comment that the name SOC is moti- 
vated by similar self-tuning to a macroscopically homo- 
geneous state under nonequilibrium conditions observed, 
for example, in "sandpile" models (Bak et ai, 1988), and 
in fluid turbulence. There, however, the self-organized 
state displays "avalanches" (vortical eddies in the fluid), 
with a power law distribution of sizes, analogous to sim- 
ilar power law-distributed fluctuations observed at equi- 
librium critical points, but without the requirement that 
an external parameter like temperature be tuned to ob- 
tain the critical state. Although the "SO" part of SOC is 
justified for the '*He state, the "C" part is not since anal- 
ogous critical power laws have yet to be demonstrated 
either theoretically or experimentally, as should become 
clear below. The name, however, has stuck and we will 
not attempt to alter convention here. 

The fact that k increases with decreasing T ensures sta- 
bility of the SOC state to small perturbations (Machta 
et ai, 1993; Ahlers and Liu, 1996). More specifically, 
an analysis of the heat diffusion equation in the normal 



**When Q and g are antiparallel (heat from below) they 
cooperate to simply produce a sharper interface. 



M 



1\ 


V. M mi 








G=0 








Q,G * 

\ G=0.04 






^ \ 
^ \ ^ 



AT 

-200 
(nK) 

-400 












-600 
-800 








100 Q 1000("W/cm' 





20 40 60 — 80 100 



FIG. 5. Simulation of the SOC state using the same simpli- 
fied 1-d model as in Fig. 4. For G oc g/Q = the temper- 
ature gradient in the normal phase [Z < 20) gives way to 
an asymptotically isothermal superfluid phase {Z > 30) with 
temperature TooiQ) < Ta.o. For G > the superfluid phase 
develops phase slips (vortices in 3-d) , and a corresponding dy- 
namic staircase structure in T{z), roughly bounded between 
Toc{z) and Tc{z), to produce the SOC state. The density 
of phase slips increases with G. Inset: experimental data 
replotted from Fig. 4 of (Moeur et al, 1997) showing the self- 
organization temperature, AT{Q) — T{Q,z) — T\{z). Only 
for Q ,$ lOOnW/cm^ (small shaded region) is the SOC state 
in the normal phase. Figure reprinted from (Weichman and 
Miller, 2000) courtesy of the Journal of Low Temperature 
Physics. 

phase (Weichman and Miller, 2000) shows that a pertur- 
bation e(x) = esoc + Se{x.) obeys an equation of motion 
whose solutions are decaying plane waves of the form 

(5e(x,t) ^e(q)e-^soc9^tg.q-(x+csoc2t)^ (4 2) 

with Dsoc = K(esoc)/C*p(esoc) and csoc = 
-|9zTA|K'(esoc)/7A,oC'p(esoc), where Cp is the equilib- 
rium specific heat at constant pressure. Thus, in ad- 
dition to the decay controlled by the diffusion constant 
i'soCi there is an unexpected anisotropic propagation 
effect where the perturbation moves upstream at speed 
Csoc- For the reasonable value Q = 50 nW/cm^ (see 
below) one finds csoc — 2.8 mm/s, and the propagation 
effect should be experimentally observable for reasonable 
cell geometries (Weichman and Miller, 2000). 

Since k diverges as e ^ 0, (4.1) impfies that esoc 
as Q — > 00. However, for large Q the Fourier law 
breaks down. In particular k{T\,Q) is finite for Q > 
(Haussmann and Dohm, 1992a), and the SOC state must 
therefore lie below Tx{z) for sufficiently large Q > Qsoc- 
This is consistent with experimental data (Moeur et at, 
1997), reproduced in the inset to Fig. 5, which shows 
that esoc < for Q > 100 nW/cm^. Using (4.1), along 
with (3.2) defining the validity of the Fourier law, one 
may show that the previous theory is valid for (Onuki, 
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1996; Weichman and Miller, 2000) 



6.1 



Q 



(l00nW/cm2 



< 1. (4.3) 



Equality in (4.3) serves an estimate for Qsoc and yields 
Qsoc — 60nW/cm2 under Earth's gravity, in very rea- 
sonable agreement with the experimental result. 

The question now remains as to the nature of the SOC 
state below Tx- The state must undergo some kind of 
transition to superfluidity, but the fact that it contin- 
ues to support a finite temperature gradient appears in- 
consistent with the isothermal nature of a superfluid. 
The resolution of this paradox is shown in Fig. 5. Let 
Too (2) = Tx{z) - Aroo(Q) and T,iz) = Tx{z) - ATe(Q) 
define local values of the interface and instability temper- 
atures discussed in Sec. III.B, where ATc = Tx^ — TdQ) 
and AToo = Txfi — Too{Q) are the deviations from Tx in 
zero gravity. For £f = an interface state, represented 
by the upper curve in Fig. 5, is formed. For g > 0, 
T(z) first drops below the local transition at a point zq; 
and attempts to asymptote to an isothermal superfluid 
at a temperature close to Toc{zo)- However, at a point 
z\ > zq, the descending line Tc{zi) meets Too (20) and 
the superfluid becomes unstable. A vortex is generated 
and crosses the cell, leading to dissipation and a finite 
temperature drop across it. By this mechanism the tem- 
perature is able to drop below the instability temperature 
and asymptote once again to an isothermal superfiuid at 
a temperature close to Toc{zi). The entire scenario then 
repeats itself approximately periodically at a sequence of 
points Zn, n = 1,2,3,..., where Tc{zn) meets Toc,{zn-i)- 
The result is a dynamic staircase structure, a snapshot 
of which is represented by the lower curve in Fig. 5. This 
structure fluctuates in time as vortices form and anni- 
hilate, but is found muiK^rically to slowly move in an 
escalator-like fashion upstream (Weichman and Miller, 
2000). 

The time resolution of present thermometry is far too 
poor to detect fluctuations in the temperature profile, 
whose mean would correspond to a straight line paral- 
lel to and somewhere between T^oiz) and Tc{z).^^ Ex- 
periments to detect the predicted stream of vortices, via 
the noise they emit in the form of second sound, are in 
the planning stages. Note that only if further numer- 
ical or experimental investigation were to reveal power 
law spatial and/or temporal correlations in the vortex 
nuclcation events (as exhibited, for example, in l//-type 
second sound noise spectra) would one obtain a postiori 
evidence in support of the "C" in SOC. 



t^A rcnormalization group calculation of this mean profile 
as a function of Q in the large-A limit is described in (Hauss- 
mann, 1999a; Haussmann, 1999b). 



A simple calculation shows that the distances be- 
tween vortices must scale as Zn — Zn-i ~ [AT^(Q) — 
ATc{Q)]ge/'yg, increasing with larger Q and smaller 
g. Although the SOC state requires finite gravity, the 
regime of widely separated vortices may in the future 
prove sufliciently interesting that a controlled low grav- 
ity experiment will become desirable. 



V. SPECIFIC HEAT AT CONSTANT HEAT CURRENT 
A. Enhanced specific heat 

The presence of a heat current is predicted to enhance 

the specific heat of superfluid ^He above its equilibrium 
value, Cq. This can easily be seen at low Q, sufficiently 
far below Tx, where the free energy enhancements are 
Ai^(r,U,) = IpsVI and A$(T,J,) = -J^/2/3,,, in 
which Ps{T) w Po|e|^ is the equilibrium superfluid den- 
sity. Thus, at fixed superfluid velocity U^, 

ACu^ = Cu^ - Co « Tp.VlCil - C)/2r^'|e|' > 0, (5.1) 

while at fixed heat current Q, 

ACq = Cq - Co « TJlCil + 0/2T^Ps\e\' > 0. (5.2) 

Closer to Tx the superfiuid density depends strongly on 
heat current and the apparent divergences at Tx in (5.1) 
and (5.2) are replaced by new singularities at the phase 
boundary Tc{Us) shown in Fig. 1(b). It is predicted that 
ACu, will remain finite, rising to a cusp at the phase 
boundary (Haussmann and Dohm, 1994), while ACq is 
predicted to diverge (Chui at al. 1996; Haussmann and 
Dohm, 1996). The latter result follows on very general 
grounds from the thermodynamic conjugacy of and 
Jg discussed in Sec. II. The usual thermodynamic ma- 
nipulations imply the relation 



Cq = Cu,+T 



fdJsV f dUs 



\dT J 



Us 



dJs 



(5.3) 



Since the susceptibility, {dUs/dJs)^ [proportional to the 
slope of the curve in Fig. 2(b)] is expected to diverge at 
the phase boundary, while {dJs/dT)us = Us{dps/dT)u, 
remains finite, Cq will exhibit a divergent enhancement. 

There have been no measurements of Cj/^ to date, but 
an experiment of this kind might be performed in the 
presence of a persistent current flowing aroTind a loop, 
similar to the superfluid gyroscope experiment (Clow and 
Reppy, 1972), where in the absence of vortices Ug indeed 
remains fixed as T is varied. A measurement coiild prove 
difficult due to the small magnitude of the enhancement 
and the challenge of holding Ug constant while measuring 
the specific heat. However, it has been suggested that 
with very fast thermometry, a measurement of Cu^ might 
be obtained through a second-sound measurement where 
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the second-sound waves are propagated perpendicular to 
XJs (Haussmann, 1997). 

Measurement of Cq is more straightforward. The pre- 
dicted divergence of ACg, together with the fact that 
Q can be experimentally controlled with great precision, 
implies a much more visible experimental signature. The 
first experimental measurements of this quantity were re- 
cently reported (Harter et al., 2000). As will be discussed 
below, they indicate that the heat capacity is indeed en- 
hanced, but with a magnitude that is significantly larger 
than theoretical predictions. 



B. The superfluid breakdown temperature 

Experiments performed at constant Q should find that 
Cq diverges at a temperature Tc{Q) < Tx, defined by 
inverting (2.20): 



kc(Q)| = 



Tx - T,{Q) 
Tx 



(5.4) 



where Qg = QoPc and x = 1/Aq = l/2i/ =i 0.746 
(Onuki, 1984; Haussmann and Dohm, 1991; Goodstein et 
al., 1996). Based on a renormalization group analysis of 
the Model F equations (Hohenberg and Halperin, 1977), 
in an approximation neglecting vortices (and hence de- 
cay of superflow), the prediction Qq w 7.4 kW/cm^ 
was obtained (Haussmann and Dohm, 1992b). More re- 
cently, using an extension of this theory, accounting for 
dissipation within a large-A^ approximation, the value 
Qq w 6.6 kW/cm^ was obtained (Haussmann, 1999b). 

These theoretical results disagree with the resiilts of 
thermal conductivity experiments (Dmican et al, 1988). 
The onset of thermal resistance was found to occur at 
a temperature which we call TbAs(Q) that obeys (5.4), 
but with X = 0.813 ± 0.012 and Qfj = 568 ± 200 W/cm^. 
This is a curve in the T-Q plane that falls below the the- 
oretically estimated TdQ) for all experimentally accessi- 
ble temperatures (see Fig. 6). A number of explanations 
have been proposed for the discrc;pancy: (a) the transi- 
tion at Tdas{Q) may be caused by a temperature insta- 
bility at the cell wall associated with the singular Kapitza 
resistance [which raises the temperature near the bottom 
(heated) end plate above that of the bulk, which there- 
fore could serve as a vortex nucleation center] , and hence 
lies below Tc{Q) (Harter et al, 2000); (b) Tdas may be 
related to a gravity-dependent transition, again lying be- 
low T'c(Q), and will increase towards TdQ) as gravity is 
reduced, e.g., by going into space (Haussmann, 1999b); 
(c) since the transition is only sharply defined in the ab- 
sence of vortices, 7bAs(Q) might be analogous to a spin- 
odal line in a first-order phase transition (Liu and Ahlers, 
1996): fluctuation-induced vortices nucleate the transi- 
tion to the dissipative phase, and Idas {Q) will diflfer from 
experiment to experiment, depending on the heating rate 



used. Since any superfluid state above Too(Q) should be 
unstable by this mechanism, an infinitessimally slow ex- 
periment should find the transition at T^{Q). 

It is possible that all of these effects (and perhaps oth- 
ers) are present. The real question, whose resolution 
clearly requires more experimental data, is which one im- 
poses the most severe limitation on present experiments. 
The answer to this question has implications for the mea- 
surement of Cq{T). If effect (c) is dominant, then exper- 
iments have basically already reached the intrinsic limit 
on how close they can approach the divergence of Cq (T) 
(for the range of Q explored thus far), though it may 
be possible to design an experiment with faster heating 
rates and fast enough thermometers to reach Tc{Q) be- 
fore a vortex can nucleate. On the other hand, if proposal 
(b) is correct, a space-based microgravity measurement of 
Cq{T) should be able to get considerably closer to Tc{Q) 
than one performed on the ground. If proposal (a) is cor- 
rect, carefully designed ground-based experiments might 
be able to approach Tc{Q) more closely: if Tdas is due to 
a boundary effect, a cell constructed with a bottom plate 
that is much larger than the cross-sectional area of the 
bulk helium sample could decrease the singular Kapitza 
resistance sufficiently so that the bulk helium can reach 
Tc{Q) without a boundary instability interfering. A cell 
of this configuration that maintains a reasonable geome- 
try for heat fiow would have to be fairly tall and might 
therefore be more susceptible to gravity effects. 



C. Experimental measurements 

The first experimental measurements of Cq(T) (Har- 
ter et al., 2000) confirm the predicted enhancement, but 
find that its magnitude is significantly larger than current 
predictions (Chui et al, 1996; Haussmann and Dohm, 
1996; Haussmann, 1999b). The data were taken over the 
range 1 /xW/cm^ < Q < 4 fiW/cm^, and a representative 
set at Q = 3.5 ^W/cm^ is shown in Fig. 7. 

The extent of the disagreement between theory and 
experiment can be observed more clearly when the data 
are plotted in scaled form. All theories predict that the 
enhancement should obey a scaling form 



ACq = |e|-«/j. 



Q 



=(r) 



(5.5) 



where fj^{x) is a universal scaling function. From (5.2) 
and the scaling relations in Table I, for a; <C 1 one has 



fM=f2X^ + 0{x^) 

/2 = [C(C + l)(QS)V2poT|5^] 



(5.6) 



Using the theoretical estimates one obtains /2 = 
8.9 J/mol K for Q3 = 7.4 kW/cm^ and /a = 7.0 J/mol K 
for Qq = 6.6 kW/cm^, where the molar volume 
27.38 cm^/mol has been used to obtain familiar units. 
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FIG. 6. Thick sohd hne: Tc{Q), the theoretically predicted 

temperature of suporfluid breakdown (Haussmann and Dolim, 
1992b). Thin solid line: rDAs(Q), whore the exponent x and 
amplitude QJj arc chosen to best fit to the observed tempera- 
ture of superfluid breakdown represented by the data points 
(Duncan et al, 1988). Dashed-dotted line: the value of Qo 
that, along with the theoretical value x = 1/2^, makes the 
experimental heat capacity data match the more recent theo- 
retical prediction for the scaling function (Haussmann, 1999b) 
(see Fig. 9 below). 

In Fig. 8 we show the heat capacity enhancement as a 
function of the scaling variable {Q/Qc)^- Since Qc{T) 
is not actually measured in the experiment, we scale 
the data using the theoretically predicted form Q^. = 
Qg|e|^« with Qg = 6.6 kW/cm^ and Ag = 2z/ = 1.342 
(Haussmann, 1999b). As anticipated, the data for all 
Q > 2 fiW/cm^ collapse onto a single linear curve, ver- 
ifying that the exponent Ag is at least consistent with 
the data. However, the slope of the experimental line is 
j^expt _ gg _j_ 4 j/jxLol K, approximately ten times larger 
than the theoretical prediction. 

An optimistic explanation for the discrepancy between 
theory and experiment is that the theories are producing 
reasonable estimates for the universal scaling function 
fj^{x), but that the nonuniversal amplitude Qq, depend- 
ing on detailed properties of the ^He system and therefore 
more difficult to compute, is estimated less accurately. 
As illustrated in Fig. 9(a), the choice = 3.4 kW/cm^W 
in fact places the experimental data on top of the more 
recent theoretical curve (Haussmann, 1999b). With this 



*''This value is in fact within the estimated margin of error 
for the amplitude calculation. The uncertainty of the theory 
is a factor < 2 (R. Haussmann, private communication), while 
the adjustment here is ~ 1.9. 



-3 -2 -1 

T-T^(^K) 

FIG. 7. Sample heat capacity data for Q = 3.5 /iW/cm^ 
(Harter et al., 2000). Thick solid line: equilibrium Co ob- 
tained from a fit to the Lambda Point Experiment (LPE) data 
(Lipa et al., 1996) that is subsequently rounded for gravity. 
Thin solid line: theoretical prediction (Haussmann, 1999b) 
that includes vortices (rounded for gravity). Solid circles: 
data from the average of the top and bottom thermometers. 
Open circles: data from the top thermometer only [beyond 
the point marked (3, where the bottom (hotter) thermome- 
ter was found to change its behavior, perhaps as a result of 
a boundary heating effect]. The temperature Tdas marks 
the onset of dissipation found in earlier thermal conductivity 
experiments (Duncan et al., 1988), while Tc(Q) is estimated 
from a certain theoretical fit to the data discussed later in the 
text. Inset: schematic diagram of the experimental cell. HRT 
stands for high resolution thermometer. 

choice, Tc{Q) lies somewhat above Tdas(Q) (vertical 
dashed line in Fig. 6). A somewhat smaller choice for 
Qq would provide an equally good fit to the earlier the- 
ory (Haussmann and Dohm, 1992b). Unfortunately, all 
of the data lie at fairly small (Q/Qc)^ < 0.3 where the 
scaling function has little structure (essentially indistin- 
guishable from linear within the scatter of the data). A 
true test would require data in the regime (Q/Qc)^ — > 1 
where the scaling function diverges. 

For completeness we also show in Fig. 9(b) an equally 
good scaling collapse based on the assumption that 
TciQ) « Tdas(Q)- Thus, we use Qc(T) derived from 
(5.4) using X = 0.813 (i.e., effectively v = 0.615) and 
Qq = 0.65 kW/cm^ [optimally chosen within the error 
bars quoted in (Dimcan et al., 1988)]. Since Tdas places 
a lower bound on Tc{Q), the sharpest conclusion we can 
make at this stage is that the data are consistent with 
the scaling hypothesis for a fairly broad range of exper- 
imentally and theoretically motivated parameter choices 
and that more data closer to Tc{Q) will be required for 
a critical test of the theory. 
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FIG. 8. Scaling plot of the differential heat capacity mea- 
surements for various values of Q (Harter et al., 2000). The 
experimental data arc scaled using Qq = 6.6 kW/cm^ and 
terminated at the temperature indicated by /3 in Fig. 7. 
The predicted collapse of the data for different Q values onto 
nearly the same curve verifies the basic scaling hypothesis 
(5.5). Thin solid line: straight line fit to the data. Thick 
solid line: theoretical prediction neglecting dissipation (Clmi 
et al, 1996; Haussmann and Dohm, 1996). Dashed line: the- 
oretical prediction including dissipation (Haussmann, 1999b). 
Neither theoretical curve is rounded for gravity. 

The experimental measurements were taken in a cell 
that was only 0.64 mm high, about as close to the op- 
timal height as practical considerations allow. Although 
the effects of gravitational rounding were therefore min- 
imized, they were not entirely eliminated. The variation 
of Tx across the cell due to gravity was ST\ ~ 8 x 10^** K. 
A reasonable criterion for a data point to be unaffected 
by gravity is that one should have \T-Tx\ > WSTx- This 
would restrict the temperature range of the experiment 
to more than about 1/xK below Tx- Essentially none of 
the interesting data shown in Figs. 7-9 satisfy this crite- 
rion. In fact, because measurements (Baddar et al., 1999) 
have shown that under a heat flux Q > 4 /xW/cm , he- 
lium exhibits appreciable dissipation, there is no range 
of parameters for which this criterion can be satisfied 
under Earth's gravity in an isothermal experiment that 
measures Cq close to Tc{Q). 

The small cell height used in the experiment precluded 
the use of side-wall thermometry. The helium tempera- 
ture was measured using thermometers mounted on the 
cell end-plates and the data were corrected for the singu- 
lar Kapitza resistance (Fu et al., 1998). The temperature 
range of the measurements was limited because the bot- 
tom (hotter) thermometer changed its behavior before 
the bulk helium temperature reached TjjAsiQ) (Harter 
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FIG. 9. Alternative scaling plot of the differential heat ca- 
pacity measurements for various values of Q (Harter et al, 
2000). The data symbols are the same as those used in 
Fig. 8. The experimental data are scaled using Qc(T) de- 
rived from (5.4) by (a) using the theoretical exponent value 
X = l/2z/ = 0.746, but amplitude Q§ = 3.4 kW/cm^ chosen 
to best match the theoretical scaling function; and by (b) as- 
suming that Tc{Q) ~ Tdas{Q) with = 0.65 W/cm^ and 
X = 0.813. For the Q-range of the experimental data, the 
two analyses are basically identical, with the lines completely 
overlapping to well within experimental resolution. Only for 
higher Q data would the two fits become distinguishable. 
Solid line: theoretical prediction neglecting dissipation (Chui 
et al., 1996; Haussmann and Dohm, 1996); Dashed line: theo- 
retical prediction including dissipation (Haussmann, 1999b). 

et al., 2000). It was proposed that this change was due 
to another boundary effect related to the Kapitza resis- 
tance. The temperature at which this phenomenon oc- 
curred is indicated by /? in Fig. 7, and is the maximum 
temperature of the data shown in Figs. 8 and 9. As a 
result of the reduced range, the data in Fig. 9 do not 
reach high enough temperatures to encounter much cur- 
vature in the scaling function. Measurements that ap- 
proach closer to the divergence are clearly needed. Data 
up to Tdas{Q) should be easily obtainable using a deeper 
cell constructed with a mid-plane thermometer. How- 
(?ver, since rounding due to gravity will be even more of 
a detriment over the region where the scaling function 
has significant curvature, the ill-effects of the deeper cell 
height will obscure some of the benefit gained by the in- 
creased temperature range. 

The only definitive way to circumvent the problems 
raisc;d in the previous two paragraphs is to perform a 
space-based measurement of Cq. An experiment in the 
absence of gravity would obtain data up to 7das(Q) 
without gravitational rounding, permitting an extension 
of the scaling data into a more revealing temperature 
range. This would allow a considerably improved es- 
timate of where Tc{Q) lies in relation to Tdas(Q) [in 
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particular if ACq is still finite at TbAs((5) one would 
conclude that Tc{Q) > TbAs(Q)]- Furthermore, a space 
experiment would test the suggestion that Tdas(Q) is 
a gravity artifact (Haussmann, 1999b), and permit the 
construction of a deep cell with asymmetric endplates to 
test whether Tdas(<3) is a boundary effect. 

VI. CONCLUSIONS 

We have reviewed the current status of experiment and 
theory of liquid helium when it is very close to the super- 
fluid transition, but not in equilibruim. Nonequilibruim 
states are most easily probed experimentally by passing a 
uniform heat flux, Q, through the liquid. If the average 
temperature of the sample is sufficiently far above Tx, 
the heat flux produces a uniform temperature gradient, 
given by (3.1). Below Tx, a sufficiently small Q produces 
a uniform temperature. Between these limits there is a 
nonlinear regime that is only partially understood. 

In the uniform temperature regime below T\, Jg (pro- 
portional to Q) and Us may be treated as conjugate 
thermodynamic variables [see (2.4)]. By analogy to equi- 
librium critical point phase transitions, one can imagine 
a critical region in the Q-T plane, and we are able to 
define new critical point exponents for these conjugate 
variables, and derive scaling laws between them that per- 
mit all of them to be evaluated in the critical region. 
These predictions have undergone only limited experi- 
mental tests. 

In all instances, the nature of the nonequilibrium phase 
transition is strongly affected by gravity. If the heat flux 
is in the opposite direction to gravity, the transition can 
be made to occur within the cell, the bottom end where 
the heat enters being normal, and the top, where the 
heat is removed, superfluid. In this case, the width of 
the interface, which is the nonlinear region, is dominated 
by Earth's gravity at the very low values of Q that are of 
primary interest to investigations of critical phenomena, 
and is too small to be probed experimentally. If the heat 
flux is in the same direction as gravity, the system devel- 
ops a temperature gradient equal to the gravity induced 
gradient in Tx, so that it has the same reduced temper- 
ature everywhere. This is the so-called SOC state. At 
small Q, the SOC state can be on the normal side of the 
transition. At larger Q it exists on the superfluid side, 
the temperature gradient being maintained by a dynamic 
mechanism of vortex production. 

On the superfluid side of the transition, the heat ca- 
pacity, Cq, is expected to be larger than Cq{Q = 0), 
and in fact to diverge at the phase boundary. Prelimi- 
nary experiments show that Cq is indeed enhanced, and 
in fact is far larger than expected theoretically. It is not 
known whether this discrepancy is due to uncertainty in 
the position of the phase boundary or to other causes. 



This result too is strongly affected by gravity, which in- 
troduces an imhomogeneity in the reduced temperature, 
e, at the small values of e necessary for these experiments. 

Although much progress has been made in understand- 
ing this nonequilibrium phase transition and exploring 
the exciting new physics that it exhibits, it is clear that 
experiments in the absence of Earth's gravity will be 
needed. One such expriment (DYNAMX) is in prepa- 
ration, and others are planned. 
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